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ABSTRACT 



^^ , Aims. We investigate the structure and shape of the photospheric and molecular layers of the atmospheres of four Mira variables. 

Methods. We obtained near-infrared A'-band spectro-interferometric observations of the Mira variables R Cnc, X Hya, W Vel, and 
RW Vel with a spectral resolution of about 1500 using the AMBER instrument at the VLTl. We obtained concurrent JHKL photometry 
using the the Mk 11 instrument at the SAAO. 
Q , Results. The Mira stars in our sample are found to have wavelength-dependent visibility values that are consistent with earlier low- 

Vh ' resolution AMBER observations of S Ori and with the predictions of dynamic model atmosphere series based on self-excited pulsation 

'tp. ' models. The corresponding wavelength-dependent uniform disk (UD) diameters show a minimum near the near-continuum bandpass 

j^ I at 2.25 yum. They then increase by up to 30% toward the HtO band at 2.0 yum and by up to 70% at the CO bandheads between 2.29 /jm 

and 2.48 //m. The dynamic model atmosphere series show a consistent wavelength-dependence, and their parameters such as the 
visual phase, effective temperature, and distances are consistent with independent estimates. The closure phases have significantly 
wavelength-dependent and non-zero values at all wavelengths indicating deviations from point symmetry. For example, the R Cnc 
^ I closure phase is 110° ± 4° in the 2.0 /jm H2O band, corresponding for instance to an additional unresolved spot contributing 3% of 

Cn , the total flux at a separation of ~ 4 mas. 

\l ■ Conclusions. Our observations are consistent with the predictions of the latest dynamic model atmosphere series based on self- 

00 ' excited pulsation models. The wavelength-dependent radius variations are interpreted as the effect of molecular layers lying above the 

^^ ' photosphere. The wavelength-dependent closure phase values are indicative of deviations from point symmetry at all wavelengths, 

|^^_ , thus a complex non-spherical stratification of the extended atmosphere. In particular, the significant deviation from point symmetry 

^— «, ■ in the H2O band is interpreted as a signature on large scales (there being a few across the stellar disk) of inhomogeneities or clumps 

in the water vapor layer. The observed inhomogeneities might possibly be caused by pulsation- and shock-induced chaotic motion in 
the extended atmosphere. 

Key words. Techniques: interferometric - Techniques: photometric - Stars: AGB and post-AGB - Stars: atmospheres - Stars: fun- 
damental parameters - Stars: mass-loss 



1. Introduction been used to characterize AGB star atmospheres. In particu- 
. , , ,.,.,, , lar, observations using the IOTA interferometer have uncovered 
Mira stars are long-period large-amplitude variable stai's on the ^^^ wavelength-dependence of Mira star diameters using a few 
asymptotic giant branch (AGB). Mass loss becomes increasingly bandpasses with spectral resolutions of up to A/AA ~ 25, pro- 
important toward the tip of the AGB, before the star evolves to ^-^^ observational evidence of molecular layers lying outside 
the planetary nebula (PN) phase, where a great diversity of mor- ^j^^ photospheric layers (e.g. ; Mennesson et al.|2002l Pen-in et al. 
phologies IS seen. Contemporary astrophysical questions related ^^^ j^^^^ molecular layers have also been present in theoret- 
to the study of Mira variables include the pursuit of the detailed ^^^j ^ -^ ^^^^^ atmospheres (Hofmann et al. mSi Ireland 
mass-loss mechanism on the AGB, including the effects of pulsa- ^^ ^j 2004a 2004b ). //-band interferometry at the IOTA interfer- 
tion and shock fronts on the structure and morphology of the ex- ^^^^^^. -^ ^^^ ^^^^^ ^^^^ (Ragland et al. 2006 ), as well as in three 
tended atmosphere and of the rtiechamsm shaping the observed ^jj^^.^ ^^^j^ ^/^^ ^ ^5 (Ragland et al. 2008 Pluzhniket al. l2009l 
PN morphologies. Near-in frared mterferometry has traditionally j^^^^ ^^^^^^^^ ^^^_^^^.^ ^j^^^^^ pj^^^^^ ^^^ ^^^^^^j ^.^^ variables, 

y, T3 J u t- J uu .u \TTn t f t r\n'n\ which reflect asymmetric brightness distributions of the photo- 

Based on observations made with the VUI Interferometer (VLTI) j b f 

at Paranal Observatory under program ID 082.D-0723 ^Phere and/or the envelope around the star. Wittkowski et al. 
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Table 1. Log of our VLTI/ AMBER observations 



Target 


Date 


JD 


Mode 


Baseline 
m 


PA 

deg 


RCnc 


2008-12-29 


4830 


MR23 


15.9/31.8/47.7 


-74 




2008-12-30 


4831 


MR21 


16.0/32.0/48.0 


-73 


XHya 


2008-12-29 


4830 


MR23 


16.0/31.9/47.9 


-71 




2008-12-30 


4831 


MR21 


15.5/30.9/46.4 


-76 


WVel 


2008-12-29 


4830 


MR23 


15.3/30.6/45.9 


-84 


RWVel 


2009-03-01 


4892 


MR23 


15.1/30.1/45.2 


-86 



Notes. The Julian date (JD) is given as JD-2450000. The AMBER 
instrument mode is MR23 (medium resolution mode with central 
wavelength 2.3 yum and range 2. 12-2.47 yum) or MR21 (2. lyum, 1.92- 
2.26 yum). The baseline is the projected baseline length for the AT VLTI 
baselines used (EO-GO/GO-HO/EO-HO) with projected baseline angle 
PA. 

Table 2. SAAO JHKL photometry 



Target 



J 
mag 



H 

mag 



K 

mag 



L Av 
mag mag 



mag 



R Cnc 0.866 
X Hya 2.375 
W Vel 2.002 



-0.155 
1.476 
0.987 



-0.652 
0.982 
0.503 



-1.063 
0.397 
-0.097 



0.07 
0.18 
0.42 



2.55 
4.12 
3.64 



Notes. The observations were obtained on 2008 December 30 (JD 
2454831), concurrently with the AMBER observations in Table[T] 



(I2008I I presented the first VLTI/ AMBER near-infrared spectro- 
interferometric observation of an AGB star, providing continu- 
ous wavelength coverage from 1.29 fiva to 2.32 //m with a spec- 
tral resolution of AjAA ~ 35. The data showed visibility and di- 
ameter variations as a function of wavelength that generally con- 
firmed the predictions of dynamic model atmospheres, where the 
diameter variations can be understood as the effects from atmo- 
spheric molecular layers (most importantly H2O, CO, and TiO). 
Here, we present the first VLTI/ AMBER observations of Mira 
variables using its K medium resolution modes with a spectral 
resolution of A/AA ~ 1500, and a comparison to the newly avail- 
able CODEX dynamic model-atmosphere series by Ireland et al. 
(i2008ii20irj . which are based on self-excited pulsation models. 

2. Observations 

We obtained /T-band spectro-interferometry with a resolution of 
A/AA ~ 1500 of the Mira variables R Cnc, X Hya, W Vel 
and RW Vel with the AMBER instrument (Petrov et al. I2007I I 
of the VLTI. Table [U shows the details of our observations. We 
obtained concuiTent near-infrared JHKL photometry of R Cnc, 
X Hya, and W Vel at the SAAO Mk II instrument. Table |2] lists 
the obtained photometry and the estimated bolometric magni- 
tudes using the procedures outlined in Whitelock et al. (2008) 
and their estimates of Ay- For R Cnc and X Hya, we obtained ad- 
ditional AMBER and SAAO observations at later epochs/phases, 
as well as coordinated mid-infrared VLTI/MIDI observations 
and VLBA observations of the SiO and H2O maser emission. 
The complete data set for these sources will be discussed in a 
dedicated forthcoming paper. 

Raw visibility and closure phase values were obtained from 
the AMBER data using version 3.0 of the amdlib data reduc- 
tion package (TatulH et al. 2007, ChelH et al. 2009). An abso- 
lute wavelength calibration was performed by correlating the 
AMBER flux spectra with a reference spectrum that includes 
the AMBER transmission curve, the telluric spectrum estimated 



with ATRAN (Lord I 19921 ). and the expected stellar spectrum us- 
ing the BS 4432 spectrum from Lan9on & Wood (2000"), which 
has a similar spectral type as our calibrators (K3-4 III). Similarly, 
a relative flux calibration of the target stars was performed using 
the calibrator stars and the BS 4432 spectrum. Calibrated visi- 
bility spectra were obtained by using an average of two trans- 
fer function measurements taken before and after each science 
target observation. The calibrators were i Hya (for R Cnc and 
X Hya), 81 Gem (for R Cnc), 31 Leo (for X Hya), GZ Vel (for 
W Vel and RW Vel), and y Pyx (for W Vel), and were selected 
from Borde et al. (I2002D using the ESO calibrator tool CalVin. 



3. Results and comparison to model atmospheres 

Figures [1] & |2] show the resulting AMBER flux, visibility, UD 
diameter, and closure phase data. The UD diameter curve is ob- 
tained by fitting a UD model to the visibility data separately for 
each spectral channel. We note that the intensity profiles are ex- 
pected to be more complex than a simple UD model, and this 
approach is used merely to readily obtain a preliminary estimate 
of the apparent diameter Signatures of a more complex model as 
well as asymmetries are most pronounced near and beyond the 
first visibility minimum, hence we used the two shorter baselines 
only for UD fits. 

The AMBER data are compared to the P/M dynamic model 
atmosphere series (Ireland et al. I2004al I2004bl) . as well as the 
newly available CODEX series (Ireland et al. l20U8l ISOTlT ). Both 
series are based on self-excited pulsation models. The only fit pa- 
rameter for each model atmosphere is the Rosseland angular di- 
ameter. Fig.[T]shows the result of the best-fit model of each of the 
P/M and CODEX series. Compared to the P/M series, the CODEX se- 
ries uses the opacity sampling method, and the model series are 
available for a wider range of basic parameters of the underly- 
ing hypothetical non-pulsating parent star The CODEX intensity 
profile was tabulated in steps of 0.0005 yum (A/AA - 40000 at 
A - 2/vm), and the synthetic visibility values averaged over the 
AMBER spectral channels assuming a Gaussian-shaped band- 
pass. The CODEX series include the series o54 (parent star pa- 
rameters: M = 1.1 Mo, L = 5400 Lo, R = 2I6R0, P = 330 
days), R52 (l.lOMo, 5200 Lo, 209 Rq, 307 d), C50 (1.35MO, 
5050 L0, 291 Ro, 430 d), and C81 (1.35 Mo, 8160 Lo, 278 Rq, 
427 d), compared to the previous series P (l.OMo, 3470 Lo, 
241 Ro, 332 d) and M (1.2 Mo, 3470 Lo, 260 Rq, 260 d). Table|3] 
provides an overview of the measured and derived quantities of 
the four targets. 

The targets of our sample have a visibility function with a 
characteristic shapqj. For each baseline, it exhibits a maximum 
near 2.25 /im and decreases towards both shorter and longer 
wavelengths. At shorter wavelengths, it shows a minimum near 
2.0/im. At longer wavelengths, between 2.3 yum and 2.5 yum, the 
visibility function displays sharp drops at the locations of the 
CO bandheads. This characteristic shape of the visibility func- 
tion of the four targets studied here is consistent with low res- 
olution (A/AA ~ 35) AMBER data of the Mira variable S Ori 
(Wittkowski et al. 2008), but shows additional details thanks 
to the increased spectral resolution. Fig.[T]displays the position 
of the H2O (1.9-2.1//m) and CO bands (2.3-2.5 fim, but where 
H2O may also contribute). 

The corresponding wavelength-dependent uniform disk 
(UD) diameters show a minimum near the near-continuum band- 
pass at 2.25 yum. They increase toward 2.0//m, where the H2O 



' The R Cnc visibility reaches into the second lobe, where the shape 
is inverted. 
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Fig. 1. R Cnc flux, squared visibility amplitude, closure phase, 
and UD diameter (from top to bottom) as a function of wave- 
length. Data of the remaining three targets are shown in Fig. 
2 of the electronic version. The wavelength range is 1.92- 
2.47 ;L/m (MR21 mode: 1.92-2.26yum, MR23 mode: range 2.12- 
2.47 jum). The modes overlap in the region 2.12-2.26/vm, where 
the two lines appear as one thicker line. The thick vertical lines 
denote the errors averaged over three wavelength intervals. The 
blue lines show a comparison to the best fitting model of the P/M 
model atmosphere series (Ireland et al. 2004a, 2004b), and the 
red lines the best fitting model of the new CODEX series (Ireland 
et al. 120081 1201 11 1. The dashed line shows a uniform disk curve 
fitted to the 2.25 yum near-continuum bandpass. 



opacity is large, by up to ~ 30%. The CO bandheads between 
2.29 yum and 2.48 jum are clearly visible for all targets with an 
increase in the UD diameter by up to ~ 70%. 

The best-fit models of the P/M and CODEX model series pre- 
dict visibility and UD curves that are consistent with the obser- 
vations. Here, the newly available models of the CODEX series 
provide a closer agreement with the data than the P/M series, at 
both the locations of the CO bandheads, which may be expected 
because of the newly introduced opacity sampling method and 
their higher spectral resolution, and in terms of the overall shape 
of the curve. The latter may be explained by the availability of 
more model series with additional parent star parameters and 
with an increased phase coverage per series. The AMBER flux 
curves are consistent with those of the dynamic model atmo- 
sphere series. In addition, the pulsation phases of the best-fit 
CODEX models are consistent with the observed visual phases, 
and their effective temperature is consistent with that derived 
from the fitted angular Rosseland-mean diameter and the bolo- 
metric flux. Furthermore, the distances obtained from the fitted 
angular diameter and the model radius are consistent with those 
derived from the measured apparent bolometric magnitude and 
the model luminosity, as well as with the period-luminosity dis- 
tance from Whitelock et al. (2008 cf. Table[3]l within 1-3 cr. 

The closure phase functions of our targets exhibit significant 
wavelength-dependent non-zero values at all wavelengths. Non- 
zero values of the closure phase are indicative of deviations from 
point symmetry, here along the projected intensity profile onto 
the orientation of the position angle given in Table[T] The consis- 
tency of the visibility amplitudes with spherical models and the 
closure phase deviation being more significant for targets that 
are well-resolved indicates that the deviation from point sym- 
metry originates from sub-structure at a relatively low flux level 
of an overall spherical intensity distribution. The strongest clo- 
sure phase signal is obtained for the most clearly resolved target 
R Cnc. It differs strongly from point-symmetry in the H2O band 
at 2.0yum with a value of 1 10° + 4°. At the continuum bandpass 
at 2.25 jum, the closure phase value is closer to with a value 
of 45° ± 5°. Towards longer wavelengths (CO band), the aver- 
age value is 153° + 4°, with values up to 165° in the bandheads, 
and thus close to a 180° symmetric value. As an illustration, the 
R Cnc closure phase value in the H2O band can be caused by 
the addition of an unresolved (up to ~ 3 mas diameter) spot to a 
circular disk, contributing up to ~ 3% of the total flux, i.e. cor- 
responding to the level of the visibility departures between data 
and spherical models, at a separation of ~ 4 mas. 



10 s 4. Discussion 



The four Mira variables of our sample exhibit consistent char- 
acteristic wavelength dependences of the visibility and conse- 
quently the corresponding uniform disk diameter that are con- 
sistent with those of earlier low-resolution AMBER data of the 
Mira variable S Ori and the predictions of the P/M and CODEX 
dynamic model atmosphere series. Here, the newly available 
CODEX series provides a closer agreement with the data than 
the earlier P/M series. This result confirms that the wavelength- 
dependent angular diameter is caused by the atmospheric molec- 
ular layers, here most importantly H2O and CO, as they are 
naturally included in the dynamic model atmosphere series. 
Concurrent JHKL photometry obtained at the SAAO was used 
to derive T^s based on the integrated bolometric flux and the fit- 
ted Rosseland-mean angular diameter. Parameters of the best-fit 
model atmospheres, such as visual phase, effective temperature. 
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Table 3. Results 



RCnc 



XHya 



WVel 



RWVel 



0V. iP) 


0.3 (356 d) 


0.7 (298 d) 


0.9 (397 d) 


(452 d) 


UD, full A [mas] (xl) 


15.2 ±0.1 (17.2) 


7.8 ±0.1 (1.7) 


10.3 ±0.1 (4.0) 


14.0 ±0.1 (6.3) 


UD, 2.25 yum (cont.) [mas] (xl) 


13.2 ±0.1 (1.3) 


6.2 ± 0.1 (0.3) 


8.0 ±0.1 (0.2) 


10.6 ±0.1 (0.3) 


UD, 2.0 //m (H2O) [mas] (^l) 


14.7 ±0.1 (4.8) 


8.3 ±0.1 (0.4) 


not obtained 


not obtained 


UD,2.4;um(CO)[mas](;r2) 


16.1 ±0.1 (8.6) 


9.0 ±0.1 (0.5) 


11.2 ±0.1 (1.0) 


14.6 ±0.1 (1.9) 


CP, full A [deg] 


110±4 


5±2 


-8 ±2 


14 ±3 


CP, 2.25 yUm (cont.) [deg] 


45 ±5 


2±2 


-4 ±2 


3±2 


CP2.0^m(H2O)[deg] 


108 ±4 


11 ±3 


not obtained 


not obtained 


CP2.4;um(CO)[deg] 


153 ±4 


5±3 


-12 ±3 


34 ±4 


Best P/M Model (xl) 


Ml In (7.4) 


P40(1.9) 


M25n (4.8) 


M12 (9.5) 


Best CODEX Model (xl) 


C50-378680 (4.6) 


O54-287940 (0.7) 


C81-250440(1.6) 


O54-262360 (4.6) 


0R„ss [mas] 


11.8 ±0.7 


4.9 ± 0.8 


5.8 ± 1.2 


9.2 ± 1.3 


yModel [j^] 


2628 ± 184 


3064 ± 380 


3216 ± 225 


3142 ±243 


'I'Modcl 


0.45 ± 0.09 


0.67 ± 0.28 


0.77 ±0.12 


0.73 ±0.10 


rM.as [K] 


2694 ± 140 


2912 ±310 


2990 ± 390 


OTboi not avail. 


dpL [pc] 


280 ± 28 


440 ± 44 


510 ±51 


366 ± 36 


^Model (0Ro.., <°*') [pc] 


252 ± 28 


465 ± 49 


362 ± 96 


217 ±33 


dModel (mbol, Mll'f') [pc] 


246 ± 46 


503 ± 80 


400 ± 37 


TOboi not avail. 



Notes. The visual phase i^vis and period P was obtained from fitting a sine function to the AAVSO and AFOEV data of the last ten periods. For 
RW Vel such data are not available, and we list the GCVS (Samus et al. 12009b period instead. The UD diameters are obtained from a fit to all data 
(full A), to the continuum bandpass 2.24-2.26 yum (2.25 yum), to the H2O bandpass 1.95-2.05 /im (2.0ytim), and to the CO bandpass 2.35-2.45 yum 
(2.4 fim). The closure phases CP were averaged over the same bandpasses. The Rosseland mean angular diameter ©ross is an average of the fitted 
diameters based on several closely fitting CODEX model atmospheres of any series with average effective temperature 7'^°'''=' and an visual phase 
•tiviodci- T^^'^" was derived from Oross and niboi (Tab.|2}. dpL is the period-luminosity distance from Whitelock et al. (2008 1. RW Vel is based on the 
same scale, but using the GCVS period and the mean K magnitude from Smith et al. J2002b . The last two rows give the distances based on Oross 



and R'^"'^"' 

Ross 



of the models at the respective phase, and based on mboi and M^^'''^'. 



and distances are consistent with independent estimates, which 
provides additional confidence in the CODEX modeling approach. 

The closure phase functions of our targets exhibit non-zero 
values at all wavelengths with a wavelength dependence that also 
correlates with the positions of the H2O and CO bands. This re- 
sult indicates a complex non-spherical stratification of the ex- 
tended atmosphere of Mira variables. The most significant devi- 
ation from point symmetry is observed in the H2O band around 
2.0yum, in particular for the clearly resolved target R Cnc. We in- 
terpret this signal as an indication of inhomogeneities or clumps 
within the water vapor layer on scales of a few resolution ele- 
ments across the stellar disk. These inhomogeneous water shells 
have also been detected for the Mira variable U Ori (Pluzhnik 
et al. I2009I I and the symbiotic Mira variable R Aqr (Ragland 
et al. I2008I I. The deviations from point symmetry at the near- 
continuum bandpass may be related either to photospheric con- 
vection cells (cf . Freytag & Hofner l2(308l l or the inhomogeneities 
of molecular layers that may also contaminate this bandpass. 

Inhomogeneous or clumpy molecular layers may be ex- 
pected as the result of chaotic motion induced by the interaction 
of pulsation and shock fronts with the extended atmosphere. The 
CODEX models indicate that the outer mass zones outward of 
1.5-2 Rosseland radii are only loosely connected to the stellar 
pulsation (cf. Fig. 1 of Ireland et al.'2008'). Likewise, outer mass 
zones on different sides of the star may be only weakly corre- 
lated with each other and may have different extensions. Icke et 
al. ( 1 1 9921) described as well that the outer layers of an evolved 
AGB star may respond with chaotic motion to the pulsations 
that originate in the stellar interior. Observations of water vapor 
layers may be particularly sensitive to these effects, but other 
molecules, such as CO, that are expected to be plentiful in the 
shocked region of the atmosphere (cf. Cherchneff 120061 ) would 



also be afl'ected by this large-scale chaotic motion. The resulting 
clumpy structure of molecular layers may explain our complex 
wavelength-dependent closure phase signal. 

Further interferometric campaigns with high spatial and 
spectral resolution are clearly needed to characterize in detail the 
morphology of atmospheric molecular layers in Mira variables. 
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Fig. 2. Flux, squared visibility amplitude, closure phase, and uniform disk diameter (from top to bottom) as a function of wavelength 
for the Mira variables X Hya, W Vel, and RW Vel (from left to right). As Fig.[T] but for the remaining targets of our sample. 



